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ABSTRACT 

In an effort to understand tlie fundamental physics of turbulent transport 
of particles and heat in a tokamak, the floating potential fluctuations in the 
the scrape-off layer plasma of ohmically heated ADITYA tokamak are anal- 
ysed for self-similarity using distribution function approach. It is observed 
that the distribution function of a sum of n data points converges to a Levy 
distribution of scale index, a = 1.111 for n < 40 and a = 2.0 for larger n. 
In both scaling ranges, the edge fluctuation is self-similar. This observation 
is backed by several supporting evidences. The results indicate that the 
small scale fluctuations transport matter and heat dominantly by convec- 
tion whereas the transport due to large scale fluctuations is by a diffusive 
process. 
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Experimental observations [1] and numerical simulations [2] have shown 
that random fluctuations in the edge plasma of a tokamak exhibit a ten- 
dency to organize themselves into coherent structures. These structures are 
expected to play an important role in causing transport of particles, a topic 
of considerable interest in modern fusion devices. The measurement of fluc- 
tuation induced particle flux shows bursty nature [3] . An attempt has been 
made to explain the bursty nature of particle transport in terms of self or- 
ganized criticality (SOC) as indicated by the observation of long range time 
correlation and self similarity in the fluctuation data [4]. However, whether 
or not the long range correlation is a manifestation of SOC (and a break- 
down of the standard transport paradigm) is an open question [5] . 

In this paper, we carry out an additional investigation of bursty nature 
of particle transport in tokamak in terms of Levy flights. The paradigm of 
Levy flight has been used recently to explain anomalous transport in such 
diverse subjects as fluid dynamics [6], transient polymers [7], subrccoil laser 
cooling [8], human heart beat [9], movement of stocks in financial markets 
[10,11] and fluid turbulence [12]. The Levy scale index, a is intimately re- 
lated to the index, v of the transport equation : < -R^(i) >= Dt^ ^ where 
< > is the mean square distance travelled by the test particle in the 

random field in time t and D is the transport coefficient. For a diffusive 
transport, v = \ whereas v > \ indicates anomalous transport. In gen- 
eral, a = 2/v and hence < a < 2 represents anomalous transport [13]. 
To demonstrate Levy process and the resultant self- similarity, we have fol- 
lowed the approach of probability distribution function (PDF) [10,11]. This 
approach is superior because the PDF involves the entire range of fluctu- 
ation amplitude unlike the rescaled range (R/S) analysis which essentially 
describes the two point function of the PDF [4]. We provide evidence of a 
self- similar non-Gaussian process, which is a Levy process, over a range of 
time scales up to 3 times the correlation time (rac)- At longer time scales the 
process is Gaussian. The observation implies a new paradigm of turbulent 
transport in tokamaks, viz., anomalous transport at short time scales and 
diffusive transport at larger time scales. 

It has been shown by Levy [14] that the sum of n independent stochastic 

variables, with a probability distribution having power-law wings, converges 
to a stable process characterized by the Levy distribution. The process 
is stable because the PDF of the sum, Z„ = Yji=iXi, of the stochastic 
variables {X}, has the same functional form as the PDF of Xj. The PDF 
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of a symmetrical Levy stable process is given by [15,16] 

1 

p{Zn) = - / exp(-n7g") cos{qZn)dq (1) 

vr Jo 

where a (0 < a < 2) and 7(> 0) are the scale index and scale factor 
respectively. If the central region of the distribution is well described by a 
Levy stable process, then the the probability of return to the origin is given 

r(i/ ) 

" ^ 7rQ;(n7)^/" ^ ^ 

where F is a Gamma function. The Levy stable symmetrical distributions 
rescale under the following transformations: 

and 

Ps{Zs) = n^/''p{Zn) (4) 
Although Levy stable process is characterized by infinite variance, the exis- 
tence of finite variance in physical systems is treated cither by introduction 
of spatiotemporal coupling [13] or by means of truncated Levy distribution 
(TLD) where an exponential tail is imposed on the values of the stochastic 
variable [15,16]. The TLD is a Levy stable process. In this paper, we re- 
port the evidence of a Levy stable process in the turbulent edge plasma of 
ohmically heated ADITYA tokamak. 

For this experiment, ADITYA tokamak is operated with with the follow- 
ing discharge parameters: plasma current, Ip = 64±4 kA, toroidal magnetic 
field at axis, Bt = 0.75 T, chord averaged plasma density, rfg = 1 x 10^^m~^, 
major radius ,Ro = 0.75 m, minor radius, a= 0.25 m. The floating potential 
is measured using a poloidal array of Langmuir probes made up of moly- 
denum wires (length = 3.5 mm, diametcr= 1 mm, separation= 3 mm) and 
located in the Scrape-off layer plasma 3 mm behind the limiter. The probes 
on a movable shaft are mounted on the top port, toroidally 72° away from 
the limiter in the electron side. The data are aquired at a sampling rate 
of 1 MHz in the time window 25-49 ms. The mean potential lies between 
-4 to 1 volts. The fluctuation data are generated after removing the 1 ms 
running average. With this procedure, we expect to remove the very low fre- 
quency fluctuations which are not part of plasma turbulence. A stationary 
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segment of data during 35-49 ms are chosen from seven discharges (8746-48 
and 8752-55) for further analysis. The autocorrelation time, Tac= 12 /xs. The 
autocorrelation function falls sharply and crosses zero at about t/Tac = 2 
with no evidence of algebric tail. 

For further analysis, stationary fluctuation time series from 7 discharges 
are normalized with their root mean square (RMS) amplitude and stacked 
together. The distribution function is obtained using 980000 data points. 
This is shown in Fig.(l), together with a Gaussian having the same vari- 
ables. The PDF of the fluctuation is nearly symmetrical (skewness = -0.38) 
but has a non-neglible kurtosis (1.24). The evidence of high kurtosis is seen 
in the distribution peaked near the origin. We next try to fit the PDF of 
the fluctuations data to a Levy distribution. 

In order to determine the Levy scale index a, the number of convolu- 
tions is chosen in the logarithmic intervals, n = 2*~^, where i =1,2, ... 11. 
The original time series is divided into non-overlapping blocks of n data 
points and a new variable, is generated for each block. The new time 
series is then subjected to Levy analysis. The probability of return to ori- 
gin, p{Zn = 0) is plotted as a function of n. This is done for the following 
reasons: (i) the distribution function is most accurate at Z„ = 0, and (ii) 
the analytical form of only p{Zn = 0) is known [Eq.(2)]. It is observed that 
on a log-log plot (Fig. 2), there are two clear slopes with values -0.9 and 
-0.5 in the n-ranges 1-40 and 40-500 respectively. The slope with value -0.9 
corresponds to Levy index a = 1.111 and that with value -0.5 to a = 2.0. 
The values of the scale factor, 7 in the two scaling ranges are 0.64±0.06 
and 8.03±0.76 respectively. It is important to note that the cross-over to a 
Gaussian process takes place at n = n^; = 40 which is much larger than that 
expected form the central limit theorem [17]. 

In Fig. (3a) we compare the PDF observed for n=l with the Levy sta- 
ble distribution of the scale index, a=l.lll and scale factor, 7=0.59. The 
Levy distribution is a good fit in the central region of the PDF upto 3a" 
values of the stochastic variable. Beyond 3a, the PDF departs from both 
Levy and Gaussian forms. Note that the values of a and 7 are obtained 
in the self-similarity range of n=l-32 ^s. This range represents small scale 
fluctuations in the floating potential. Figure 3(b) shows the full distribution 
function for the n— values 1, 2, 4, 8, 16 and 32. As the number of con- 
volutions increases, the distribution function becomes broad and the peak 



4 



value (or, the probability of return to origin) decreases. If the stochastic 
variable, Z„ and the distribution function, p{Zn) are rescaled in accordance 
with Eq.(3) and Eq.(4) respectively, the distribution functions, PsiZs) col- 
lapse on n =1 distribution (Fig. 3(c)). Thus, the distribution functions 
of different convolutions are self-similar in the range n =1 to 32. A self- 
similarity is also obtained for n— values 64, 128, 256 and 512 by using a=2.0 
and 7=8.03ib0.76. Figure 4 shows the rescaled distribution functions in the 
two scaling regimes on the same graph. The difference of log{ps{Zs = 0)) 
values in the first and the second scaling regimes provides a measure of the 
'distance', A between the two scaling regimes [15]. The A ^ 0.6 for the 
results shown in Fig. 4 indicates that the 'distance' is significantly large. 

We have verified the above results by carrying out the following anal- 
ysis: (i) the root mean square (RMS) amplitude of the stochastic variable 
{Zn), az„ shows a power-law behaviour, az„ ~ n*^/^. The exponent, z/=1.8 
in the short time scaling range (n=l-32) whereas, iy=l.O in the long time 
scaling range (n= 100-800). Thus, the small scale fluctuations show non- 
Gaussian PDF and a self-similar scaling behaviour. The long scales show 
self-similarity of Gaussian type. It can be argued that the RMS amplitude 
of potential fluctuation is proportional to the root mean square distance 
travelled by a test particle in the random field [18]. (ii) the rescaled range 
{R/S) analysis on the combined data set shows a somewhat longer scaling 
range at short time scales (n= 1-128) having Hurst parameter, H=0.9. At 
long time scales, n=256-4096, the value of H = 0.6. One is tempted to 
interpret it as an evidence of long range self-smilarity of non-Gaussian type. 
However, it should be pointed out that R/S range analysis has tendency to 
give somewhat higher values of H [19]. Therefore, we consider that at long 
time range, the fluctuation data show self-similarity of Gaussian type, (iii) 
when we calculate the kurtosis (K) of the time series corresponding to dif- 
ferent convolutions (n) , it is observed that kurtosis decreases with increasing 
n and K < 0.2 for n > 60. This result is similar to those reported earlier 
in tokamak [20] and stcllerator [21]. Thus, the small scale fluctuations are 
non-Gaussian whereas the large scales are Gaussian. 

In conclusion, we have presented evidences of a Levy stable process in 
tokamak edge turbulence. The sum of n data points of a fluctuation time 
series converges to a Levy distribution of scale index, a= 1.111 for n < 40 
and a= 2.0 for larger n. The probability distribution functions are self sim- 
ilar in both scaling ranges. This observation is backed by several supporting 
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evidences. 

We finally give a possible interpretation of these observations and spec- 
ulate on their implication for the problem of bursty transport due to turbu- 
lence in the tokamak edge region. We have earlier carried out conditional 
statistical analysis of Langmuir probe data in the edge of tokamak ADITYA 
(for similar discharges) and shown that the fluctuations are dominated by 
coherent structures which appear intermittently in space and time, typically 
lasting 25-30 (^s and having poloidal scale lengths of a few cm. These struc- 
tures are also associated with sharp radial potential gradients separating 
the last closed magnetic surface from the scrape-off layer plasma. Bursts 
of turbulent transport into the scrape-off layers are related to intermittent 
breakdown of radial confinement and the coherent structures extending into 
the open field line regions. We believe that the non-gaussian PDFs of po- 
tential fluctuation (^) described above are due to such coherent structures. 
For time scales longer than the life time (~ 25-30 /is), we observe only a 
randomized average behaviour and hence infer only Gaussian statistics and 
diffusive behaviour. For shorter time scales, on the other hand, convective 
effects due to coherent structures dominate and we observe the anomalous 
behaviour, < >~ with i^=1.8. We may also make an estimate of the 
enhancement of turbulent transport due to convective effects introduced by 
the presence of coherent eddy-like structures by using the analysis of Rosen- 
bluth et al. [22]. The enhancement factor / = D*/D ~ P^/^ where the 
Peclet number P ~ vd/D is the ratio of the diffusion time to the eddy turn- 
over time. For a typical eddy size of d ~ 1 cm, the typical eddy velocity, 
V ~ cE/B 10^ cm/s and diffusion coefficient, D ~ 10^ cm^/s [1,23], we 
find an enhancement factor / ~ 3. Thus, the presence of transient coherent 
eddies can give sudden avalanche- like enhancements of transport by a factor 
of 3 or more, resulting in bursty nature of transport. 
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FIGURE CAPTION 



(1) The probability distribution function (PDF) of floating potential fluc- 
tuation. The amplitude (Z) is normalized to the standard deviation. 
The symbol (open circle) shows the experimental data points whereas 
the dotted line shows a Gaussian PDF of unity standard deviation. 

(2) The probability of return to the origin, p{Zn = 0) as a function of the 

number of data points (n) in the non-overlapping blocks. For n=1024 
and beyond, the PDF is not well defined and p{Zn = 0) estimates may 
be in error. Since we have filtered out time scales larger that 1 ms, 
they do not contribute to the PDF. 

(3) (a) The PDF, p{Z) for the experimeental data (bullets) together with 
the Gaussian distribution (dotted line) and the Levy distribution (solid 
line), (b) the PDF, p{Zn) vs. Z„ for different n. The width of p{Zn) 
increases with increasing n, (c) the rescaled PDF, Ps{Zs) for the first 
six convolutions (n =1-32). 

(4) The comparision of rescaled PDF, Ps{Zs) in the scaling ranges, n =1- 
32 (solid lines, without symbol) and n =64-512 (with symbols). 
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Figure 1 
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Figure 3 
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Figure 4 
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